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ABSTRACT 

We investigate the formation of the stellar halos of four simulated disk galaxies using high resolution, 
cosmological SPH + N-Body simulations. These simulations include a self-consistent treatment of all 
the major physical processes involved in galaxy formation. The simulated galaxies presented here each 
have a total mass of ~ 10 12 M Q , but span a range of merger histories. These simulations allow us to 
study the competing importance of in-situ star formation (stars formed in the primary galaxy) and 
accretion of stars from subhalos in the building of stellar halos in a ACDM universe. All four simulated 
galaxies are surrounded by a stellar halo, whose inner regions (r< 20 kpc) contain both accreted stars, 
and an in-situ stellar population. The outer regions of the galaxies' halos were assembled through 
pure accretion and disruption of satellites. Most of the in-situ halo stars formed at high redshift out 
of smoothly accreted cold gas in the inner 1 kpc of the galaxies' potential wells, possibly as part of 
their primordial disks. These stars were displaced from their central locations into the halos through 
a succession of major mergers. We find that the two galaxies with recently quiescent merger histories 
have a higher fraction of in-situ stars (~ 20 — 50%) in their inner halos than the two galaxies with 
many recent mergers (~ 5 — 10% in-situ fraction). Observational studies concentrating on stellar 
populations in the inner halo of the Milky Way will be the most affected by the presence of in-situ 
stars with halo kinematics, as we find that their existence in the inner few tens of kpc is a generic 
feature of galaxy formation. 

Subject headings: Galaxy: formation — Galaxy: halos — galaxies: formation - galaxies: halos - 
methods: N-Body Simulations 



1. INTRODUCTION 

Accreted stellar halos are a natural consequence of 
galaxy formation in a ACDM universe, in which the 
smallest structures collapse first an d galaxy formation is 
thought to proceed b ottom-up (e.g. IWhite fc Rees|[l978t 
iSearle fc Zinnl Il978| ). Halos, composed of both dark 
and baryonic matter, grow by merging with other ha- 
los. While the gas from mergers and accretions loses 
its energy through cooling and settles into a disk, the 
accreted stars and dark matter form a halo around the 
galaxy. Galactic stellar halos are thus predicted to be 
built from multiple accretion events starting from the 
first structures to collapse in the Universe. 

It is for this reason that much work has focused on 
using kinematically identified halo stars to look for fossil 
records of the galaxy's accretion history. The possibility 
remains, however, that stars formed within the potential 
well of a galaxy can become displaced from the inner- 
most regions into a kinematic stellar halo. In order to 
draw conclusions about the accretion history of a galaxy 
based on halo stars, one must know how abundant such 
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in-situ stars are in the halo, and how they may dilute the 
signatures of accretions and mergers. 

Owing to the low surface brightness of halos, the only 
presently well studied stellar halos are those of the Milky 
Way and M31. Current observational studies of the 
Milky Way's halo find strong evidence that its outer 
regions (r > 20kpc) were at least partially assembled 
through accretion events. Numerous stellar streams be- 
lieved to have originated from small accreted galaxies 
have been observed i n the Milky Way' s outer halo (e.g. 
Newberg et al.ll2002t lYannv et al.ll2003t iBelokurov et all 
20061 : iGrillmairl l2009h . Multiple wraps of stars believed 
to have been tidal l y removed from the Sa gittarius dSph 
(jlbata et all 119941: iMaiewski etal] I2003D highlight the 
process of disrupting satellites in the Milky Way's po- 
tent ial, and serve as s trong evidence of hierarchical merg- 
ing. iBell et all (|2008( ) have found, using F and G turn off 
stars in SDSS, that the distribution of halo stars is highly 
structured, and the outer halo is consistent with being 
formed almost entirely from accreted subhalos. M31's 
stellar halo is observed to contain possibly even more 
substructure and st reams than what is observed around 
the M ilky Way (e.g. llbata et al J200lL [2007t iGilbert et all 
2007), suggesting its recent merger history may have 
been more active. Despite the plethora of streams ob- 
served around the Milky Way and M31, both the fraction 
of halo that is accreted, as well as the number of distinct 
accretion events leading to the observed substructure, is 
unclear. 

Although the nearby halo of our Galaxy has been stud- 
ied more extensively than the outer halo, in some ways 
its properties are less well understood. This is due in 
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part to the short dynamical timescales close to the disk, 
where accreted stars would no longer betray their origin 
through spatial structures. It is in these inner several 
tens of kpc of the halo that one might expect dissipa- 
tive processes to contribute to the overall stellar popu- 
lation. Observational evidence of the contri bution of ac- 
cretio n to the inner halo is slowly growing. iHelmi et al.l 
(1999) have detected two streams in the solar neigh- 
borhoo d by studying the ang ular momentum of stellar 
orbits. iMorrison et alj (|2009f) have also recently shown 
that the angular momentum distribution of a sample of 
stars with a median height above the Galactic plane of 
2 kpc does not seem to be smooth, possibly indicative 
of an accretion origin. A hybrid theory in which the 
Milky Way was assembled through a combination of hi- 
erarchical accretions and dissipativc processes has been 
used by som e to u nderstand prop erties of field halo stars 
e.g. lNorrislH99l iCarollo et all l200l . iChiba fe Beersi 
2000), using a large set of halo stars selected with no 
kinematic bias, find a vertical gradient in the mean ro- 
tation of metal poor s tars, a signature of dissipational 
processes. Ilbata et al.l (|2007l) . using number count maps, 
have shown that a large fraction of the halo of M31 ap- 
pears to be rather smooth. The authors interpret these 
observations to suggest that some of the halo of M31 may 
have formed through dissipative processes. However, the 
degree to which dissipation is important in the forma- 
tion of stellar halos, or one of their components, remains 
uncertain. 

Simulations have also been extensively used to in- 
ves tigate the formation of ga lactic halos. The work 
of iBullock fc Johnston! (|2005f ) combines dark matter 
only N-body simulations with semi-analytic models to 
study the build up of stellar halos via accretions alone. 
Their work has been highly successful in explaining the 
discrepant chemical abundances o f Milky Way sat el- 
lite galaxies and fie ld halo stars ([Font et al.l l2006al lbl: 
Uohnston et al.ll2"008| ). While this powerful and efficient 
tool is particularly suited to better understand progen- 
itors and signatures of accretion events, it can not tell 
us about the possible relative importance of dissipative 
proces ses in the assembly of a galactic halo. iBrook et al.l 
(|2003f) have used N-Body + SPH to simulate the for- 
mation of disk galaxies in a self consistent way, follow- 
ing the evolution of dark matter, gas and star forma- 
tion. They find that recently accreted satellites con- 
tribute stars to the halo which can be distinguished from 
an early formed h alo through phase-space information. 
lAbadi et~a l. ( 2006), using a set of fully cosmological sim- 
ulations, find that in-situ stars dominate the inner 20 kpc 
of their galaxies, while accreted stars make up the ma- 
jority of the outer regions. While the Abadi work distin- 
guishes between the accreted and in-situ populations in 
their galaxies, the work does not investigate the possible 
importance of in-situ stars in simulated stellar halos. 

In this paper, we use four high resolution cosmologi- 
cal SPH + N-Body simulations to investigate the origin 
of stellar halos. These self-consistent simulations permit 
the study of the competing importance of in-situ star 
formation and accretion from subhalos in the formation 
of stellar halos in a ACDM universe. Furthermore, the 
supernova feedback and cosmic UV background imple- 
mented in our simulated galaxies results in a luminous 
satellite population that is more similar to the Milky 



TABLE 1 

PROPERTIES OF SIMULATED GALAXIES 



Run 


M vir 

Mq 


Ntot at z=0 
dark+gas+stars 


particle 

Mq 


M* 

particle 

Mq 


kpc 


MWlhr 


1.1 x 10 12 


4.9 X 10 6 


7.6 x 10 5 


2.7 x 10 4 


0.3 


MWlmed 


1.1 x 10 12 


6.1 x 10 5 


6.1 x 10 6 


2.2 x 10 5 


0.3 


Gall 


3.3 x 10 12 


3.7 x 10 6 


2.6 x 10 6 


9.2 x 10 4 


0.3 


b.277 


7.4 x 10 11 


2.3 x 10 6 


1.2 x 10 6 


4.6 x 10 4 


0.35 


h285 


7.7 x 10 11 


3.0 x 10 6 


1.2 x 10 6 


4.6 x 10 4 


0.35 



a Gravitational softening length 

Way's than those in previous N-Body + SPH studies. 
This improvement is essential to model realistic stellar 
halos. We quantify the degree to which the halos are ac- 
creted, and note the effect of accretion histories on the 
overall formation of the stellar halos. While past numer- 
ical efforts have concentrated on the growth of stellar 
halos through pure accretion, we show that both accre- 
tion and in-situ star formation contribute to the inner 
regions of stellar halos. 

The paper is organized as follows. Section 2 describes 
the details of the simulations used, and the methods used 
to determine the origin of halo stars. In Section 3 we dis- 
cuss the dual origin of the stellar halos of our simulated 
galaxies, and investigate how the properties of the in-situ 
and accreted populations in the stellar halos vary with 
merging history. Section 4 discusses possible numerical 
and resolution effects on the formation of the stellar ha- 
los in our simulations, as well as the effect of feedback on 
the accreted fraction. Section 5 discusses the connection 
to observational results and concludes. 

2. SIMULATIONS 

The high resolution simulations used in this study 
were run to z = 0, as part of an ongoing simulation 
project. The two more m assive galaxies (MWlhr, Gall) 
are described in detail in iGovernato et alj ([2007, here- 
after G07), but have since been rerun at eight times the 
mass resolution and twice the spatial resolution. The 
reader is referred to G07 for the full details of these sim- 
ulations. The remaining two lower mass galaxies (H277, 
H285) are a continuation of t hese efforts. H277 is de- 
scribed in iBrooks et al.l (|2009l ). We briefly summarize 
the salient features of these simulations here. 

The halos used in this study were selected, based 
on their mass and merger history, from low resolu- 
tion, dark matter on ly simulations run with GASOLINE 
(fWadslev et al.ll2004D . using WMAP year 1 cosmology 
for MWlhr and Gall, and WMAP year 3 cosmology 
for H277 and H285. Each halo was then resimulated at 
higher resolution using the v olume renormalization tech- 
nique (|Katz fc Wh ite 1993). This approach simulates 
only the region within a few virial radii at the highest 
resolution, while still simulating a large enough box to 
account for the effect the large scale tidal field has on the 
angular momentum of the halo. 

We concentrate a large part of this work on MWlhr, 
the highest resolution N-body + SPH simulation of a 
Milky Way mass galaxy run to z — published to date. 
We single this simulation out for extensive study as it has 
a lower resolution counterpart, run with the identical ini- 
tial conditions, but 1/8 the mass resolution of MWlhr, 
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as well as several runs completed with different treat- 
ments of feedback. The previous extensive study done on 
MWlhr (see G07) permit the numerical and resolution 
tests necessary for this work. MWlhr has 4.9 million 
particles (dark matter + gas + stars) within its virial 
radius, and a mass of 1.14 x 10 12 M Q . We define the 
rcdshift of last major merger as the time at which a sec- 
ondary galaxy, with a mass no less than 1/3 of MWlhr's 
mass, first enters the virial radius of MWlhr. This oc- 
curs at a redshift of 4, although the merger is not com- 
plete until z ~ 2.5. The merger history of MWlhr is 
thus similar to that of the Milky Way, which is thought 
to have experienced its last major merger at z ~ 2.5 
(jHammer et al.ll2007f ). The mass of star particles in this 
simulated run is ~ 2.7 x 10 4 M Q , high enough to re- 
solve an in- falling dwarf galaxy of stellar mass 10 9 Mq 
with more than 10 4 star particles. A physically moti- 
vated recipe was us ed to describe star f ormation and su- 
pernovae feedback ([Stinson et al.l 120061 ) . with a uniform 
UV backgro und that turns on at z=9 , mimicking cosmic 
reionization (jHaardt fc Madaul fl996). The properties of 
all the simulated galaxies used in this work are described 
in Table 1. 

All dark matter halos and subhalos in the simula- 
tion are identified using AMIGA'S Halo Finder 7 (AHF, 
iKnebe et al.l l200"lh . Each halo is identified using a min- 
imum of 64 dark matter particle mem bers, above whic h 
the halo mass function converges fG07. lReed eTaII[2003l) . 
The properties of each simulation particle are output at 
z=6, z=5, z=4, and every 320 Myr after until z=0, for a 
total of 43 time steps output. 

2.1. Selection of Halo Stars 

Before we can study the mechanisms that built-up the 
stellar halos of our simulated galaxies, it is first neces- 
sary to identify the star particles which make up the 
galactic halos. In both our study, as well as in observa- 
tional work, it is important that the stellar samples used 
to draw conclusions about halo properties and forma- 
tion are not contaminated by a population of disk stars. 
Unless an observational sample only includes stars that 
reside several disk scale heights above a galaxy's mid- 
plane, the use of kinematic and metallicity information 
is necessary to obtain a clean sample of halo stars. Some 
observational works rely on samples of stars chosen based 
on a metallicity cut, where the calculated [Fe/H] is sev- 
eral tenths of dex below that of disk s tars (for example 
iBeers et al.ll200d : IChiba k Beersll2000H . Other observa- 
tional studies have used the proper motions of stars t o 
isolate those with halo kinematics (S mith et al.l l2009f) , 
while o ther works combine both metallcities an d kine- 
matics (|Carorlo et al.ll2007t iMorrison et al.ll2009ft . 

We use the full simulation information to identify the 
most complete and uncontaminated set of halo stars pos- 
sible. We have chosen to define our halo sample using the 
full 3 dimensional phase space information available for 
all stars in the simulations. Rather than implementing 
a specific spatial cut, this approach allows us to study 
the detailed origin of all stars belonging to a halo. If we 
had instead chosen to define the stellar halo of each sim- 
ulated galaxy as the population of stars with distances 

7 AMIGA Halo Finder, ava ilable for download at 
|http://www. aip.de/P eople/aknebe /AMIGA] 



greater than ~ 2 scale heights above the galactic disk, 
the trends we report in the following sections would be 
unchanged. However, our results then would not have 
been applicable to the many detailed observational stud- 
ies of the Milky Way's halo that are focused on stars 
within d < bkpe. In a future paper, we will investigate 
the possible observational consequences of this predicted 
origin. 

A kinematic analysis of the star particles in each galaxy 
was done to decompose them into disk, bulge, and halo 
kinematic components. To decompose the disk from the 
spheroid, we first align the angular momentum vector of 
the disk with the z-axis. This is done in order to place 
the disk in the x-y plane. We then calculate the angular 
momentum of each star in the x-y plane, J 2 , as well as the 
momentum of the co-rotating circular orbit with similar 
orbital energy, J C irc- In this definition, stars with circu- 
lar orbits in the plane of the disk will have J z / Jdrc ~ 1- 
Disk stars are selected with a cut of J z /J C irc > 0.8, which 
is equivalent to an eccentricity cut of e < 0.2. This con- 
dition matches t he eccentricities obsery ed of disk stars in 
the Milky Way (|Nordstrom et al.ll2004h . 

Once the disk population has been identified, we de- 
fine a sphere centered on each galaxy, with a radius equal 
to that of the disk's. For stars within this inner sphere, 
the spheroid (halo + bulge) of each galaxy is defined us- 
ing a cut in J z /J C irc such that the net rotation of the 
spheroidal population equals zero. The spheroidal com- 
ponent of the galaxy is further decomposed into halo and 
bulge because a break in the mass density profile is ob- 
served. Stars whose total energy is calculated to be low, 
and hence are tightly bound to the galaxy, are classi- 
fied as part of the bulge. All other stars are placed in 
the halo category. The energy cut between the halo and 
bulge is set in order to make a two component fit of 
the mass profile with bulge stars within the "break" and 
halo stars outside the "break" . Bulge contamination is 
not a concern because the energy cut also ensures that 
stars which reside close to the center of the galaxy, and 
are classified as belonging to the halo, have orbits which 
take them far from a classical bulge. Moreover, obser- 
vational studies that rely on kinematics to select halo 
stars would have associated our halo stars with a halo, 
not bulge, component. We show in Section 3 that the 
stellar populations we associate with the halo extend out 
to regions far beyond what would be considered a bulge. 
Stars which reside beyond the inner sphere defined by 
the disk's radius are all classified as halo stars, regard- 
less of their kinematics. We have verified that the halo 
population has a rotational velocity with a Gaussian dis- 
tribution. 

The Jz/Jdrc cuts described above leave a set of stars 
whose kinematics do not make either the disk or spheroid 
cuts. These stars are dominated by thick disk and pseudo 
bulge stars. We note that using kinematic information 
for decomposing galactic components unavoidably leads 
to some slight overlap in populations. 

For the remainder of the paper, we use the term "halo 
stars" to refer to stars identified using this kinematic 
definition. There is no bias for our selected halo in-situ 
stars to be prograde in any galaxy. This indicates that 
we have not mistakenly classified a significant number of 
disk stars (in particular thick disk stars) as halo stars in 
our analysis. 
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2.2. Halo Star Origin 

We trace each star particle located within the main 
galaxy at z = back in time. We identify at each time 
step the dark matter (DM) halo to which each particle 
belonged using AHF. Stars are considered bound to a DM 
halo if they are identified as being a member of that halo 
for two consecutive timesteps. We follow the particles' 
histories back to z = 6, as well as follow the gas particles 
from which stars have formed. This procedure results 
in four different classifications for stars' origin: accreted, 
in-situ, ambiguous, and other. These classifications are 
described in the following sections. 

2.2.1. Accreted vs. Insitu 

We use these stellar histories to distinguish between 
stars that were accreted and those that formed "in-situ" . 
Accreted stars formed inside of halos other than the pri- 
mary dark matter halo, but through accretion and merg- 
ing have become unbound to their progenitor and now 
belong to the primary halo. In-situ stars formed from 
gas particles bound to the primary. These stars formed 
within the primary potential well and will be referred 
to as the "in-situ" sample. In Section 3.1 we show that 
while these in-situ stars reside in the halo at z=0, they 
formed near the center of the galaxies' potential wells, as 
part of their primordial disks. The gas particle progen- 
itors of these in-situ stars ended up in the primary halo 
from both smooth g as accretion, particularly in cold flows 
(|Brooks et al.| [2009i , as well as from subhalos that under- 
went gas stripping as they merged with the primary. All 
galaxies are traced back far enough in redshift, to z = 6, 
that we can accurately determine the star formation his- 
tory of the vast majority (> 98%) of halo stars. Approxi- 
mately 1% of the stars in these simulations formed before 
z=6, and as such we can not determine their formation 
origin. We can thus accurately distinguish the accreted 
stars from the in-situ stars for the vast majority of par- 
ticles in each simulated run. 

2.2.2. Ambiguous Origin Due to Time Resolution 

Our simulations also contain stars whose formation ori- 
gin is unclear due to the limited number of time steps 
output for each simulation. For example, if a gas parti- 
cle spawns a star particle at approximately the same time 
that it first becomes bound to the primary, it is uncertain 
whether that star formed in the primary or in its origi- 
nal subhalo. Such stars are referred to as "ambiguous" 
in the reminder of the paper. These stars comprise less 
than 10% of all the kinematically identified halo stars in 
all of the simulations (see Table 2). 

2.2.3. "Other" Stars in the Simulation 

To ensure that the in-situ stars are a robust feature of 
these simulations, we now investigate whether this pop- 
ulation of halo stars were formed preferentially near the 
density and/or temperature thresholds for star forma- 
tion, where spurious star formation may occur due to 
numerical effects. The spatial resolution of these simula- 
tions is ~ 300 pc, much coarser than the scales on which 
star formation occurs in real galaxies. It is therefore 
necessary to use a prescription for star formation that 
captures globally the observed trends in galaxies. The 
prescription used here for star formation enforces a lo- 
cal gas density greater than 0.1 cm -3 and a temperature 



colder than 30, 000 K at which stars form. 8 If a gas parti- 
cle matches these criteria, there is some probability that 
it will form a star particle. This probability is scaled 
to match the observed Kennicutt -Schmidt law for star 
formation (see lStinson et aLll2006l for further details). 

We found that in-situ stars that formed at a look back 
time of 10 Gyrs or more did not form near the den- 
sity and temperature cutoffs. They formed under the 
same conditions as stars with similar ages in the disk 
of the galaxy. However, in-situ halo stars in the simu- 
lation younger than 10 Gyr old formed at preferentially 
lower densities and higher temperatures than their young 
disk counterparts, near the implemented cutoff for den- 
sity threshold. While we have no reason to believe these 
late forming in-situ stars are a numerical artifact, chang- 
ing our prescription for star formation would possibly af- 
fect whether or no t these in-situ star s had formed at all 
in the simulation. iMaller fc Bullock! (|2004h have shown 
that accounting for hot halo gas instabilities may be im- 
portant in understanding the total baryonic mass of a 
galaxy. These authors show that the fragmentation of 
the hot halo gas results in hot low density gas which 
does not cool, and is able to survive in the halo of a 
galaxy. This might explain the young in-situ stars in our 
simulation that formed from a low density hot gas. To be 
maximally conservative on our prediction of the extent 
to which in-situ stars compose the halo, we have not in- 
cluded in-situ stars that formed within the last 10 Gyr in 
any of the following analysis. The broad conclusions of 
this paper are not affected by whether or not we choose 
to ignore these intermediate and young in-situ stars. We 
refer to these stars as "other" for the reminder of the pa- 
per, and only consider insitu stars formed more than 10 
Gyr ago when discussing "in-situ" stars. The total mass 
contributed by other and in-situ stars to each simulated 
galaxy is listed in Table 2. 

3. DUAL ORIGIN FOR STELLAR HALOS: A GENERIC 
FEATURE OF GALAXIES 

The presence of both in-situ and accreted stars is a 
generic feature of the kinematically defined halos of all 
four simulated galaxies. Each simulated stellar halo was 
found to be composed of a purely accreted outer halo, 
and an inner halo where both accreted and in-situ pop- 
ulations reside. This is illustrated by Figure [T] The red 
solid lines in this figure show the fractional contribution 
of in-situ stars to the stellar halos as a function of dis- 
tance. The blue dotted lines and green dash-dot lines 
show the same for stars accreted to the halo at look- 
back times of more and less than 9 Gyr, respectively. 
This figure shows that all four stellar halos have an in- 
situ component that primarily resides within their inner 
~ 20 kpc, and that the accreted component of each stel- 
lar halo extends to more than 100 kpc from the center of 
each galaxy. The outermost regions of the stellar halos 
of all four simulated galaxies are dominated by stars ac- 
creted within the last 9 Gyr. The central location of all 
four in-situ populations shows that their presence is most 
relevant for observational studies that are limited to the 

8 This temperature floor is due to our adopted cooling curve, 
which assumes primordial abundances without metals, does not 
yet include molecular cooling, and thus cannot extrapolate below 
~10,000 K. It is assumed that dense gas at this temperature will 
continue to rapidly cool. 
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inner halo of the Milky Way. Table 2 lists the fractional 
contribution, as well as the total mass, of the in-situ and 
accreted populations in each of the four simulated halos 
in our sample. 

3.1. In-Situ Star Formation 

We go on to study in detail the formation properties 
of in-situ stars in the halo of MWlhr. As explained in 
Section 2, we focus on this particular simulated halo be- 
cause of the extensive numerical and resolution testing 
which it has undergone. We find that the in-situ stars of 
Gall, h277, and h285 formed and populated their stellar 
halos by the same qualitative process described below for 
MWlhr. 

In order to understand where the in-situ stars in the 
halo of MWlhr originated, we studied the gas particles 
that formed this stellar population. The vast majority of 
the gas progenitors were brought into the primary galaxy 
in smooth gas flows. Fewer than 2% of the in-situ stars 
were formed from gas stripped from accreted subhalos. 
Cold gas flow was found to also b e the dominant cont rib- 
utor to the disk stars of MWlhr ([Brooks et alJl2009f ). 

To investigate how in-situ stars came to form part of 
the stellar halo of MWlhr, we first study the locations 
at which these stars formed relative to the center of the 
galaxy. Figure 2 shows the radial distribution of in-situ 
stars at the time of their formation relative to the center 
of the primary (black solid line). This shows that in-situ 
stars form near the very center of MWlhr's potential 
well. 

Approximately 70% of the in-situ stars were formed by 
z ~ 3. Between 2 < z < 3, MWlhr experienced at least 
3 significant mergers with mass ratios ^{ MWlhr < 10, de- 
fined at the time of virial radius crossing. Substantial 
mergers disrupt the orbits of stars in both the primary 
and merging galaxy. As two galaxies of approximately 
equal mass collide, their potentials change rapidly, caus- 
ing the energies of their stars to change as well. The 
violent relaxation process will cause some stars to gain 
energy, and alter their initial orbits. We find that by z=2 
(red dotted line in Figure 2), the distribution of in-situ 
stars is no longer sharply peaked at the inner 1 kpc of the 
galaxy, but is instead peaked at 5 kpc from the galactic 
center. This transformation from centrally located orbits 
to halo orbits was most likely due to the major mergers 
which MWlhr experienced between 2 < z < 3. After 
z=2, MWlhr does not undergo any more mergers with 
such small mass ratios, and so we expect that the ra- 
dial distribution of in-situ stars will not change. Indeed, 
we find that their distribution at the present day (blue 
dashed line in Figure 2) is very similar to that at z=2. 

The origin of the halo in-situ stars in cold gas flows, 
and their formation in the innermost regions of MWlhr's 
potential well all highlight that this population of stars 
formed in a rapid inflow and collapse of gas that also 
formed the disk of MWlhr, but were displaced through 
merger events into the halo. The process by which these 
insitu stars formed is reflected in that their relative im- 
portance in each stellar halo rises toward the center of 
each simulated galaxy. 

3.2. Comparing The Accreted Halos With Other Works 

Much work has been done recently with dark mat- 
ter only N-body simulations coupled with semi-analytic 



models to study the growth of Milky Way-like galactic 
halos via pure accretion of stars. In this section, we com- 
pare results from such numerical studies with our findings 
for MWlhr's accreted halo. We focus here on MWlhr 
as it was chosen to have a total mass, and merger his- 
tory similar to that of the Milky Way ([Governato et al.l 
l2007h . making it a good candidate to compare with pre- 
vious works which have also focused on simulating Milky 
Way- like stellar halos. We find that the accreted halo of 
MWlhr displays properties which are typical in compar- 
ison with other numerical studies. 

As shown in the top left panel of Figure 1, accreted 
stars dominate the halo of MWlhr at all radii. While 
the innermost 20 kpc of the halo contains both accreted 
and in-situ stars, the halo beyond this is highly domi- 
nated by accretion. The earliest accreted stars, which 
became bound to the primary halo of MWlhr at a look- 
back time of ~ 9 Gyr or more, dominate the accreted 
population of the halo out to 30 kpc from the center. In 
the halo's outer regions, however, the majority of stars 
were accreted more recently. 

The accretion history of MWlhr 's stellar hal o agree s 
well with the numerical study of iHelmi et al.l (|2003h . 
These authors found that ~ 90% of stars are in place 
in the inner 10 kpc of a Milky Way-like halo 10 Gyr 
ago, while the mass grow th of the outer reg ions is more 
continuous. The work of iFont et al.l (|2006af ) also found 
that nearly 80% of the stars in the inner 20 kpc of sim- 
ulated halos were in place at a lookback time of 9 Gyr. 
I Font et aTl (|2006bl ) have shown in their numerical work 
that the inner regions of halos are typically assembled 
from a few satellites with stellar masses M* ~ 10 9 M Q . 
In MWlhr, the majority (75% by mass in stars) of the 
accreted mass in the inner 50 kpc originated in two sub- 
halos, each with total mass ~ 2.5 x 10 10 M Q . 

The accreted stars which reside in the halos of Gall 
and H285 differ significantly from those found in the stel- 
lar halo of MWlhr. Unlike the accreted population of 
MWlhr, the majority of the accreted stars in the ha- 
los of Gall and H285 became bound to their respective 
stellar halo more recently, as shown in Figure 1. This is 
likely due to the fact that both Gall and H285 have had 
a more active recent merger history, as described in the 
following section. The accreted stellar halo of H277 was 
assembled at similar times to that of MWlhr, although 
the accreted population of this stellar halo contributes 
far less to the overall mass of the halo in comparison to 
MWlhr. 

3.3. Possible Role of Merger History 

Despite the universal presence of both in-situ and ac- 
creted halo stars in all of the simulations, the relative 
contributions of these dual components varies widely 
among the four stellar halos. The fractional contribu- 
tions of these stellar populations to the simulated stellar 
halos, as well as their total masses, are summarized in 
Table 2. The fraction of accreted halo stars ranges from 
~ 30 — 85%, while the fraction of in-situ halo stars ranges 
from less than 10% to more than 50%. The high and low 
ends of these ranges do not sum to 100% because of the 
small presence of the ambigious stars. Although the frac- 
tional contribution of in-situ stars covers a wide range, 
the total mass of these stars in all four simulations is 
substantial. The mass of in-situ stars in both the stellar 
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Fig. 1. — The relative contribution of accreted and in-situ stars to the stellar halos of the simulated galaxies. Accreted stars are separated 
by the time at which the stars first became bound to the primary halo. The blue dotted line, and the green dash-dot line show the accreted 
stars which became bound to the primary galaxy more than 9 Gyr ago, and less than 9 Gyr ago, respectively, while the red solid line shows 
the in-situ stars. Stars whose origins are unknown, as described in Section 2.2, were omitted in this figure, and their fractional contribution 
ignored. 



halos of H277 and Gall is ~ 4 x 10 9 M Q , and - 1 x 10 9 M Q 
for MWlhr and H285. 

All four of these simulated galaxies have similar total 
masses, however they span a range of merger histories. 
Three of the galaxies, MWlhr, H277, and H285 each host 
a stellar disk with a mass of ~ 2 x 10 10 Mq, while the 
mass of Gall's stellar disk is ~ 6 x 10 10 M Q . The total 
halo mass of H277 and MWlhr is 6 - 9 x 10 9 M Q , and 
- 3 x 10 10 M Q for Gall and H285. The bulge to disk ra- 
tio (B/D) of each simulated galaxy was obtained using a 
two component Sersic + exponential fit to the one dimen- 
sional radial profile of the I band surface br ightness maps 
made with the SUNRISE software package ( Jo nsson et all 
[20061) . The B/D for MWlhr, Gall, H277, and H285 are 
0.37, 0.66, 0.29, and 1.2, respectively (Brooks et al, in 
prep). What role might merger history play in the range 
of accreted/in-situ fractions observed in these four sim- 
ulations? There is not an obvious trend between epoch 



of last major merger and either in-situ fraction or aver- 
age assembly time of the accreted component; Gall and 
MW1 had their last 3:1 merger at z ~ 4 and H277 and 
H285 at z ~ 2.5. However, there appears to be a possible 
trend with the galaxies' recent merger histories. Figure[T] 
illustrates one manifestation of this. The accreted stel- 
lar halo components of MW1 and H277, the two galaxies 
with the highest in-situ fraction, were largely assembled 
more than 9 Gyr ago, with the innermost regions of these 
halos containing very few stars accreted in the last 9 Gyr. 
However, Gall and H285, which contain few in-situ halo 
stars, have accreted stellar halos that, on average, were 
assembled much more recently at all radii, with at least 
30% accreted within the last 9 Gyr. It is important to 
recall that we have defined accretion time as the time 
at which a star becomes unbound from its subhalo, and 
bound to the primary dark matter halo, which is a dif- 
ferent, and invariably later, time than when the subhalo 
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TABLE 2 
STELLAR HALO ORIGINS 



Run total halo: accreted: fraction (%) in-situ: fraction (%) ambiguous: fraction 3 "other" stars: b 

mass (Mq) mass (Mq) mass (Mq) mass(M0) mass (Mq) 



MWlhr 




70 


21 


9 






6.22 x 10 9 


3.8 x 10 9 


1.14 x 10 9 


4.99 x 10 8 


7.8 x 10 8 


MWlmed 




66 


25 


9 






5.24 x 10 9 


3.23 x 10 9 


1.21 x 10 9 


4.42 x 10 8 


3.57 x 10 8 


MWlthermal c 




85 


12 


3 






9.74 x 10 9 


7.96 x 10 9 


1.19 x 10 9 


2.17 x 10 8 


3.6 x 10 8 


Gall 




85 


12 


3 






3.39 x 10 10 


2.84 x 10 10 


4.17 x 10 9 


1.06 x 10 9 


4.04 x 10 8 


H277 




32 


58 


10 






9.03 x 10 9 


2.7 x 10 9 


4.85 x 10 9 


7.89 x 10 8 


6.9 x 10 8 


H285 




87 


7 


6 






2.93 x 10 10 


2.2 x 10 10 


1.71 x 10 9 


1.56 x 10 9 


3.9 x 10 9 



a The origin of these stars is ambiguous due to the finite time resolution of the simulation, as described in Section 2.2.2 
k These stars are excluded from the analysis of this paper, as discussed in Section 2.2.3 

c MWlthcrmal was run with the same mass resolution as MWlmed. but with a different supernova feedback recipe. This is 
discussed in further detail in Section 4.1 
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Fig. 2. — The radial distribution of MWlhr's in-situ halo stars 
at the time of their formation (black solid line), at z=2 (red dotted 
line), and at present day (blue dashed line), relative to the center 
of MWlhr. 

entered the primary's virial radius. 

In order to more explicitly study the merger histo- 
ries of the simulated galaxies, we traced the total num- 
ber of luminous satellites that enter the virial radius 
at different redshifts. In particular, we are interested 
in those satellites whose masses are in the range 10 < 
M, 



primar 



V/M: 



satellite < 200 at the time of merging, as 
these are expected to be the buildi ng blocks of the ac- 
crete d component of stellar halos (|Bullock fc Johnston! 
2005). Table 3 shows the number of such merged satel- 
lites at redshifts smaller than z = 2. The stellar halos 
of Gall and H285, the two galaxies with more recently 
assembled stellar halos, have continued to accrete a large 
number of satellites after z=2. The masses of these 
galaxies' accreted stellar halos are comparatively large 
(> 10 10 M Q ), overwhelming the population of halo in-situ 
stars. In contrast, MWlhr and H277, whose stellar halos 
were found to have assembled at more ancient times, have 
had a more quiet recent merger history. These galax- 



ies merged with fewer satellites in the studied mass and 
rcdshift range, and have the less massive accreted ha- 
los. MWlhr and H277 also both have halos which host 
a comparatively larger in-situ population than Gall and 
H285. 

4. NUMERICAL AND RESOLUTION STUDIES 

Owing to a dual origin, the detailed properties of a 
stellar halo are a function of both i) the properties of the 
primary galaxy itself, and ii) the properties of the pro- 
togalaxies that merged to form its accreted component. 
Accurate modeling of a galaxy's stellar halo thus requires 
an accurate treatment of the formation and evolution of 
that galaxy over a wide dynamical range. In this section 
we discuss the effects of feedback, as well as numerical 
and resolution effects, on the overall formation of stellar 
halos in our simulated galaxies. 

4.1. The Effect of Feedback on the Simulated Galaxies 

The four simulated galaxies used in this study have 
sizes, ages, and metallicities that match tho s e ob - 
served at low redshift. G07 and iGovernato et"ai1 ()2008l ) 
showed that disk galaxies naturally form in a cosmo- 
logical context in simulations such as those studied 
here, and have sizes that place them on the I-band 
Tully-Fisher re l ation (TF e.g. iGiovanelli et a .1 119971 



iMcGaughl [20051 : iGeha et al.ll200l . iBrooks et all (f2007h 
used a set of GASOLINE simulations (including MW1) 
to show that galaxies with total z = halo masses in 
the range 3.4 x 10 9 M© < M* < 1.1 x 1O 12 M have 
(O/H) abundances th at match both those measured by 
iTremonti et~aT1 12004) for more than 53,000 SDSS galax- 
ies loc ally, and those measu red by Erb et al. (2006) at 
z — 2. iMaiolino et a l. (2008) demonstrated that galaxies 
simulated with GASOLINE (including the progenitors of 
H277 and H285) are also a good match to the observed 
mass-metallicity relationship at z = 3. 

The success of these simulations in matching the ob- 
served mass-metallicity relationship at va rying redshifts 
(jBrooks et al.1 120071 : IMaiolino et aill2008l ) demonstrates 
that these simulations have made substantial progress 
in overcoming the historic "overcooling problem" (e.g. 
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TABLE 3 
LOW z MERGERS IN 
SIMULATED GALAXIES 



Run mass ratio mass ratio 
[10:50] [50:200] 



MWlhr 2 5 

Gall 2 10 

H277 1 2 

H285 7 8 



Note. — The number of merged 
satellites within the given mass 
ratio for each simulated galaxy 
(M pr i marB /M sateiHtc ), at times 
more recent than a redshift of 2. 



iMaver et alj [2008) . Gas in simulations has tradition- 
ally cooled rapidly, forming stars quickly and early, and 
thus producing a mass mctallicity relation that is too 
enriched at a given stellar mass, particularly at high z. 
In the simulations examined here, supernova energy is 
deposited into the ISM mimicking the blast-wave phase 
of a su pernova (following t he Sedov- Taylor solution, de- 
tails in lStinson et al.|[2006f ). Cooling is turned off for gas 
particles within the supernova blastwave radius for a pe- 
riod of time after the supernova event. Supernova feed- 
back regulates star formation efficiency as a function of 
halo mass, resulting in the mass-metallicity relationship 
described above. This regulation of star formation also 
leads to realistic trends in gas fractions, with our low- 
est m ass galaxies being the most gas rich (jBrooks et al.l 
l2007f) . and reproducing the observed incidence rate and 
column densities of D amped Lyman a systems at z = 3 
(jPontzen et al.ll2008f l. 

Feedback, with its role in regulating star formation, 
will also effect the luminosity function of galaxies in a 
simulation, which will in turn effect the properties of 
the accreted stellar halo. For example, if the simulated 
dwarfs have too many stars relative to dwarf galaxies in 
the real universe, then we may overestimate the relative 
importance of accreted versus in-situ stars in the halos. 
We therefore compare the luminosity function of satel- 
lites within the virial radius of our two medium resolu- 
tion galaxies, MWlmed and MWlthermal. We use these 
two simulated runs because they have the same mass 
resolution, but different implemented feedback mecha- 
nisms. MWlmed adopts the blastwave model discussed 
above, and used for all of the high resolution runs studied 
here. In MWlthermal, supernova energy was deposited 
as thermal energy to the surrounding gas particles but 
cooling is not turned off in these particles. These gas 
particles will then quickly radiate away the deposited en- 
ergy, and this simulation will suffer from the overcooling 
problem. 

Using the age and metallicity of each stellar parti- 
cle in a satellite, we calculated the absolute V band 
magnitudes of the z = satellites of MWlmed and 
MWltherma l using the stellar population models of 
Starburst99 (jLeitherer et al.lfl~999t I Vazquez fc Leithererl 
2005). Figure 3 shows the luminosity function of the 
satellites in MWlmed and MWlthermal, as well as the 
observed lum i nosity function of Milky Way dwarfs, from 
iGrebel et ail {2003). The satellites of MWlmed are sim- 
ilar in total number and luminosity to the Milky Way's, 



except at the bright end of the luminosity function. The 
satellites of MWlthermal, however, are too numerous 
at all luminosities. The discrepancy between the lumi- 
nosity function of MWlmed and MWlthermal is due 
to the overcooling problem typical of simulations with 
thermal feedback, which causes satellites to form which 
are highly centrally concentrated and contain too many 
stars. Because the MWlthermal run produces over- 
luminous dwarf galaxies, we expect that its halo will 
be more dominated by accreted stars than the halo of 
MWlmed. As shown by Table 2, which lists the accreted 
and in-situ fractions of both of these galaxies, the halo 
of the MWlthermal has ~ 20% more accreted stars in 
its halo relative to in-situ stars than MWlmed (85% vs 
66%) and < 1/2 of the fraction of in-situ stars (12% vs 
25%). While the total mass of in-situ halo stars is about 
the same for both simulated runs (1.2 x 1O 9 M0), the total 
mass of accreted stars in the thermal run is 8 x 10 9 M©, 
and only 3.2 x 10 9 M Q in the run with blast-wave feed- 
back. 

In summary, several features of the simulations we an- 
alyze here contribute to bringing the number of lumi- 
nous satellite galaxies into good agreement with those ob- 
served around the MW: Our uniform UV background un- 
binds baryons from simulated halos with total masses be- 
low ~ 10 9 M Q , creating "dark" halos and bringing down 
the number of luminous satellite galaxies. Our supernova 
feedback scheme helps to significantly reduce the num- 
ber of stars produced in the remaining low mass satellite 
galaxies. However, because our satellites are still slightly 
too bright (ie, have too many stars), our accreted frac- 
tions are an upper limit. This only serves to increase 
the potential relative contribution of in-situ stars in the 
inner halos of galaxies. 

4.2. Resolution Effects on Star Formation 

iBrooks et al.1 (|2007l ) found that several thousand par- 
ticles were required inside of the virial radius of a halo 
in order for the star formation history to converge. For 
the high resolution simulations studied here, this conver- 
gence criterion is met in halos with total masses above a 
few xl0 9 M Q . For the medium resolution runs, the con- 
vergence limit is met in halos with total masses above 
2 x 10 10 M©. Resolution testing has shown that poorly 
resolved halos have smaller total stellar masses than ha- 
los which are fully resolved. We therefore underesti- 
mate the contribution to the primary stellar halos from 
dwarf galaxies with total mass below this convergence 
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Fig. 3. — The luminosity function of satellites for two simulated 
runs with different feedback prescriptions, as well as the observed 
lumin osity function of the Milky Way's dwarfs, from fGrebel et all 
p003h . 

limit. However, we find that dwarf galaxies less mas- 
sive than 10 10 M Q contribute less than 5% to the overall 
mass in accreted stars in the high resolution stellar ha- 
los. Due to the resolution effects discussed above, this 
small fraction of accreted satellites will contribute even 
fewer stars at lower resolution than they would if they 
were fully resolved, decreasing their overall contribution 
to the accreted stellar halo of the medium resolution runs 
(MWImed and MWlthermal). This effect is minimal 
due to the tiny contributio n of these halos. Our find - 
ings agree with the work of Bullock & Johnston (2005), 
who showed that satellites with total mass greater than 
2 x 10 10 M Q contribute 75 - 90% of the mass to the pri- 
mary stellar halo. 

Because the in-situ stars are forming deep in the po- 
tential well of the high resolution main halo, we do not 
expect their formation properties to vary with resolution 
because the star formation history of the main halo is 
fully resolved. In order to test this assumption, we com- 
pared the in-situ formation in the high resolution MWlhr 
run with the lower resolution MWImed run. We expect 
that if the formation of these stars was caused by a nu- 
merical resolution effect, then their fractional contribu- 
tion to the simulated halos would change as the resolu- 
tion at which the galaxy is simulated changed. We find, 
however, that at both resolutions, in-situ stars make up 
a substantial fraction of the stellar halo of the Milky Way 
massed galaxies. Table 2 shows that the overall in-situ 
fraction of the stellar halos of MWlhr and MWImed are 
21% and 25%, repectively. The 4% larger contribution 
of in-situ stars to the stellar halo of MWImed is likely 
due to the slight underestimating of accreted stars in this 
medium resolution run, as described above. 

5. DISCUSSION & CONCLUSIONS 

In this paper we have analyzed four SPH + N-Body 
simulations of approximately L* disk galaxies to investi- 
gate the formation of their stellar halos. We have shown 
that the inner halos of all four galaxies contain both stars 
accreted from merged satellites, as well as in-situ stars 
formed within the primary galaxy. While these results 
are theoretical, they are backed by recent observational 
studies of the stellar halos of the Milky Way and M31. 

Observational work has shown that the properties of 
Milky Way halo s tars exhibit signatu res of a possible 
double population. [Car olio et al.l (|2007f ) have found that 



a sample of local halo stars in SDSS can be thought of as 
two distinct, but overlapping, components, with different 
spatial distributions, orbits and metallicities. They at- 
tribute, as others have, these different stellar properties 
to different formation mechanisms for the two halo com- 
ponents. They conclude that the inner halo may have 
formed dissipatively from the gas rich merger of early 
sub-galactic halos, while the outer halo stars were ac- 
creted from subhalos which merged with the Milky Way's 
dark matter halo. 

Observations of M31 have revealed that this giant spi- 
ral galaxy is also surrounded by a large extended stellar 
population. M31 has been shown to have a distinct stel- 
lar halo with a surface brightness profile similar to the 
Milky Way's, and a metallicity gradient with the inner 
regions of its h alo (R< 30 kpc) more metal rich than the 



outer regions (jGuhathakurta et~aT1 120051: iKalirai et aTl 
l2006f) . Using number count maps. llbata et al!l ~i 2007f ) 



find a large metal poor smooth component to the halo 
of M31, with a more metal rich inner halo. These works 
highlight the possible signatures that M31's stellar halo 
may, like the Milky Way, have two underlying compo- 
nents, with different formation histories. 

The results of such observational findings for the ha- 
los of the Milky Way and M31, however, have previously 
lacked counterparts in simulations as most of the nu- 
merical studies done to date on stellar halos have con- 
centrated on their build up through pure accretions. Our 
work begins to bridge the gap by studying self-consistent 
simulations where the contribution of in-situ stars to the 
overall halo have been taken into account. Our main 
results are: 

• In-situ stars are a generic feature of kinematically 
defined stellar halos. These stars formed deep in 
the potential well of the primary galaxy and were 
later displaced to halo orbits as a result of mergers. 

• In-situ stars reside primarily in the inner few tens 
of kpc of the galaxies' halos. 

• Stars accreted off satellites dominate the stellar 
halo at all radii for three of the four galaxies stud- 
ied. The accreted stellar halo components were 
found to extend to out to a few hundred kpc. 

• The fractional contribution of in-situ stars to the 
simulated galaxies' halos varies greatly, from 5 — 
50%. 

• The two galaxies in our sample whose halos had the 
smallest contributions from in-situ stars were also 
those with accreted halos that assembled less than 
9 Gyr ago. This is likely because such galaxies have 
more active recent accretion histories, which lead 
to massive accreted halos. The large contribution 
of recently accreted stars overwhelm the underly- 
ing in-situ populations in these halos. The other 
two galaxies in our sample, with large contribu- 
tions from in-situ stars, had a significant fraction 
of their accreted stars in place at times earlier than 
9 Gyr ago. The relative importance of in-situ stars 
in these halos is due to the recent quiet merger his- 
tory of these galaxies, which have led to less mas- 
sive accreted halos. 
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• Different supernova feedbacks greatly impact the 
formation of simulated stellar halos. This under- 
scores the necessity of carefully implementing feed- 
back to produce a realistic satellite luminosity func- 
tion, and to hence conduct realistic studies of halo 
formation. 

• The overall importance of in-situ stars to the stellar 
halos found in this study is likely a lower limit. This 
is because i) we have ignored a fraction of in-situ 
stars for numerical reasons, and ii) the satellites 
of the simulated galaxies studied here are brighter 
than those observed in the Milky Way. The calcu- 
lated contributions of accreted stars to the stellar 
halos are thus an upper-limit. 

• The qualitative trends of our results are unchanged 
if we had instead used a spatially defined set of 
stellar halos. 

Our results suggest that galaxies with more active re- 
cent mergers, like M31, may host halos where the relative 
importance of an in-situ population is low. Galaxies, 
such as the Milky Way with little evidence for recent 
merger events, on the other hand, may have a larger rel- 
ative contribution from an in-situ population in its inner 
regions, along with its accreted stars. However, a larger 
set of simulations with a wide range in merger histories 
is needed to fully understand the complex relationship 
between halo formation and mergers. 

Vast observational resources are currently dedicated to 
the study of the Milky Way's and M31's halos. Spatial, 
kinematic and chemical abundance information of halo 
stars being collected and used to attempt to reconstruct 
the accretion and formation history of the galaxies. Our 
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results highlight that such data sets, particularly ones 
focused on local, and hence inner, halo stars, will con- 
tain a significant population of stars born in-situ in the 
galaxies. These in-situ stars may dilute the signatures 
of accretion history in observations. Given the short dy- 
namical timescales in the inner halo, one may not expect 
to be able to discriminate these two populations using 
their kinematics alone. However, given the very different 
sites of their formation, the chemical abundances of in- 
situ and accreted stars may be quite distinct. In a future 
paper, we will investigate possible ways to observation- 
ally disentangle the two halo populations. 
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